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Abstract. We consider a Lorentz-violating theory of gravity where the aether vector is 
taken to be nondynamical. This “ponderable aether theory” is almost the same as Einstein- 
aether theory (where the aether vector is dynamical), but involves additional integration 
constants arising due to the loss of initial value constraints. One of these produces an effective 
energy density for the aether fluid, similar to the appearance of dark matter in projectable 
Hofava gravity and the mimetic dark matter theory. Here we investigate the extent to 
which this energy density can reproduce the phenomenology of dark matter. Although it is 
indistinguishable from cold dark matter in homogeneous, isotropic cosmology, it encounters 
phenomenological problems in both spherically symmetric configurations and cosmological 
perturbations. Furthermore, inflationary considerations lead us to expect a tiny value for 
the ponderable aether energy density today unless a sourcing effect is added to the theory. 
The theory then effectively reduces to dynamical Einstein-aether theory, rendering moot the 
question of whether an aether must be dynamical in order to be consistent. 
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1 Introduction 

Some recent attempts to address problems posed by cosmology and quantum gravity involve 
imbuing spacetime with additional structure. When such structure violates local Lorentz 
boost symmetry, it singles out a preferred time direction at each point, defining an aether 
field In a companion paper [1], we considered a variety of ways in which the aether can 
be constructed from background structures, i.e. fields that appear in the action but are fixed 
rather than varied in the action principle. Here, we analyze some of the phenomenological 
consequences of these theories in astrophysical and cosmological settings. 

We will focus on two particular structures from ref. [1], the fixed threading and the 
fixed clock. A threading is a congruence of curves that foliate a manifold [2, 3]. In the 
fixed threading theory, the aether lies tangent to a threading that is fixed as background 
structure. This fully determines the direction of the aether vector. A clock is a scalar field 
Tp that determines the overall scale of the aether vector via a'P = 1- When the aether 
is constrained to be a unit vector, ip corresponds to proper time along the flow lines of 
hence its interpretation as a “clock.” 

The aether theories constructed using these background structures resemble Einstein- 
aether theory, but with the crucial difference that some or all components of tt“ are non- 
dynamical, i.e. some of the aether equations of motion are not explicitly imposed. In the 
“least dynamical” situation, with both a hxed threading and a fixed clock, none of the aether 
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equations of motion are imposed, and the dynamics is solely determined through the Einstein 
equation. Nevertheless, the Bianchi identity guarantees that the Lie derivative along of the 
Einstein-aether field equation vanishes. The theory is therefore equivalent to Einstein-aether 
theory with a constant source in the aether equations of motion. 


The present work looks into the phenomenology of these sources in the aether field 
equation. The most notable effects come from ponderable aethers, which arise in fixed clock 
theories. The corresponding source, fl, functions as an energy density in the Einstein equa¬ 
tion. In cosmological settings, it contributes a dust-like fluid component to the Friedmann 
equation. This raises the question of whether it can account for the observed cold dark 
matter (CDM). We investigate this idea on both small scales, involving the Newtonian limit 
and static spherically symmetric solutions, as well as large scales, involving linear cosmo¬ 
logical perturbations. Both regimes possess significant phenomenological difficulties. There 
are no static, spherically symmetric solutions with nonzero values for the source densities, 
precluding a steady final state for spherically symmetric gravitational collapse. What hap¬ 
pens instead remains an open question. For cosmological perturbations, the source density p, 
leads to rapidly growing scalar perturbations on superhorizon scales, which is at odds with 
observations of cosmological structure. 


In Ref. [1], we noted that the ponderable aether theory is closely related other theories 
where dark matter arises as an integration variable. Notably, when the aether is constrained 
at the level of the action to be geodesic and hypersurface-orthogonal, the result is the IR 
limit of projectable Hofava gravity [4]. The “dark matter as an integration constant” [5] in 
that theory is analogous to the source density p for the ponderable aether theory. Another 
similar theory with such a dark matter component is the mimetic matter theory [6], which 
was shown to be equivalent to Einstein gravity together with a constraint that the gradient 
Va<p of a scalar held be a unit one-form [7]. The arguments of Refs. [8, 9] then show 
that the mimetic theory is equivalent to projectable Hofava gravity with the GR choice of 
coupling parameters. Since projectable Hofava gravity arises as the limit c^j —>■ oo, Ca —>■ oo 
of Einstein-aether theory [10, 11] {c^j, Ca are dehned in (2.1b)), the results of this paper apply 
to the mimetic and projectable Hofava theories by taking this limit. Note also that the vector 
mimetic theory of Ref. [12] is closely related to the ponderable aether, although the aether 
is dynamical in the former but nondynamical in the latter. 


This work begins in section 2 with an overview of the theories described in [1], and 
discusses the dynamics when one or more of the aether equations are not imposed. In partic¬ 
ular, the constraints for the theory, described in 2.2, require special attention. In section 2.3, 
we show how the ponderable aether theory arises as a limit of Lorentz-violating CDM [13]. 
Following this, we proceed to phenomenology on small scales in section 3, where the obstacle 
coming from static, spherical symmetry is discussed. Finally, we explore cosmological solu¬ 
tions in section 4, focusing particularly on the challenges raised by the scalar perturbations. 
Section 5 provides a summary and discussion of the results. 


Throughout this paper we employ the metric signature (-I-, —, —, —). Latin letters denote 
abstract tensor indices, while Greek letters are used for coordinate indices. Tildes are used 
to represent tensor densities of weight -|-1. 
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2 Lorentz-violating structures 
2.1 Dynamics for Lorentz-violation 

The coupling between and gab is described through an action that incorporates all terms 
consistent with the unbroken spacetime symmetries. When organized in a derivative ex¬ 
pansion, the lowest order terms appear at two derivatives, and comprise the Einstein-aether 
action [14, 15], 

Sa, = J d'^xy/^ {R + K), (2.1a) 

K = ^9“^ + CaO-"^ + + CaO^. (2.1b) 

The terms in (2.1b) are the expansion 9 = VaxR-, acceleration at = u^'VaUb, shear Uab = 
V \9fab - U(^aO‘b) {fab = UaUb “ gab is the Spatial metric) and twist iOab = V [aUb] - U[„ab] 
of the aether congruence. The coupling constants appearing here are related to the usual 
constants ci,2,3,4 by = ci -|- 3c2 + C3, c^- = ci -|- C3, = ci - C3, and Ca = ci -|- C4 [11]. In 

addition, the unit constraint on can be imposed by adding a Lagrange multiplier term 

S\ = - j (fxyf^X {gabu'^u^ - 1 ) . ( 2 . 2 ) 

Finally, the coupling to matter helds y is described by a matter action Sm[Xi9ab\ which 
involves no direct interaction between ri“ and y. The total action for the theory is then 

S = S^ + Sx + S^. (2.3) 

The various constructions of described in [1] have the effect of producing theories 
where different components of can be considered nondynamical, i.e. parts of its equation of 
motion are not imposed. Nevertheless, when the Einstein equation and matter field equations 
hold, the aether field equations are almost imposed; specifically, their derivative along the 
flow of vanishes, 

£u {Sa - 2\Ua) = 0. (2.4) 

s s 

Here the covector density £a = ^ is the functional derivative of the aether action (2.1a) 

with respect to u°‘ (see equation (A.l) for its explicit form), and A = y/—gX. This relation 

follows from the diffeomorphism invariance of the action (2.3), and is equivalent to the 
conservation of the aether stress tensor. 

These nondynamical aether theories can therefore be thought of as Einstein-aether the¬ 
ory supplemented with a source for the aether equation of motion, 

£a - 2A'Ua = -jla, (2.5) 

where the source must be invariant along the flow of u°‘, 

£u£a = 0 . ( 2 . 6 ) 

For different types of background structures, different components of jla are nonzero. In the 
fixed threading theory, where the aether is constructed as tangent to a nondynamical threading 
of spacetime by timelike curves, the quantity fia = i^a ~ u^xLa)fi,b is nonzero. Alternatively, 
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in the fixed clock theory, the parallel component of the aether is determined by a clock field 
background structure by requiring u ■ dfi = 1. In this theory, the quantity 

h = (2.7) 


is nonzero. The fixed aether theory possesses both a fixed threading and a fixed clock, which 
together completely determine u“. All components of {la are generically nonzero for the fixed 
aether theory. 

The Einstein equation obtained from varying with respect to the metric takes the form 

^ab + ^UaUb = — —Tab, ( 2 - 8 ) 


1 ss 

where T„h is the matter stress tensor, and T®? = , 

ab ^gab 

holding hxed (see equation (A.2)). Equation (2.5) can 
fj., in which case (2.8) becomes 


is the metric variation of (2.1a) 
36 used to eliminate A in favor of 


^ab 2^ ^c^aUb + ^flUatlb — (2.9) 

The left hand side contains the usual Einstein-aether metric equation terms plus an additional 
energy source coming from /r. The aether internal energy is conserved, in the sense that 
is a conserved current density: 


Va(/ltt“) = £u{l = Tuiu^'fla) = 0. 


( 2 . 10 ) 


Thus {1 is like the density of a conserved mass that is carried along by the aether flow. 
The presence of this additional mass density that follows flows of motivates the name 
ponderable aether for these theories. 


2.2 Constraints 

Einstein-aether theory is a diffeomorphism invariant theory, which implies that some of its 
equations of motion are constraints rather than dynamical equations. These constraints place 
restrictions on the allowed initial values in the theory. 

The constraints involve a combination of the metric (2.8) and aether (2.5) equations of 
motion [16, 17]. Eor a given time coordinate t (a scalar function on spacetime) the constrained 
quantities in the absence of matter are 

= Vat + u'^Sb) . (2.11) 

This combination of the equations of motion involves one fewer t-derivative than the other 
equations, hence their role as constraints on initial values. 

Since these quantities involve the aether held equation, they will not vanish in the 
nondynamical aether theory. Instead, the source density fib characterizes the constraint 
violation, 

Cf) = -{£ut)ilb. (2.12) 

The evolution of the constraint violation in time is determined by (2.6), which together with 
(2.12) implies 

£uCP = -i£a£ut)yib. (2.13) 
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For generic t this is not a simple evolution equation for however, if t is a proper time 
coordinate r for i.e. = 1, then the right hand side of (2.13) vanishes. In an 

aether-adapted coordinate system, where the aether has components 


u“ = (l,0,0,0), (2.14) 

the constraint quantities satisfy 

= Oi (2.15) 

i.e. they are simply constant along the flow of u°‘. 

Another situation where the time evolution of the constraint quantities simplifies is 
when t = r] is the parameter for a vector field proportional to the aether, with 

N a scalar function and = 1- This example is relevant in cosmological applications, 

in which r] is the conformal time and N = a is the scale factor. In this case the constraint 
evolution is given by 


= -{£^N-^)£b - N-^£^£b 
= -{N£uN-^ + N-^£uN)£h - £Jb - 

= —jl'S/bliiN. (2-16) 


Contracting this equation with yields the evolution of the scalar constraint, 

N£^{uXi^^) = 0 . 


(2.17) 


In a coordinate system where C 
an overall 77 -derivative, 


drj is the time evolution vector, this equation reduces to 



(2.18) 


and allows for immediate integration. Furthermore, if N depends only on 77 and not on 
the spatial coordinates, the spatial components of (2.16) vanish. This means the spatial 
constraints also satisfy the simple evolution law, 


^ Mn) 

dr] * 


= 0 . 


(2.19) 


When we examine cosmological perturbations in section 4.2, we choose a gauge where 
the aether is aligned with the time evolution vector dr^ and receives no linear perturbations. 
In this case, N = a{r]) is independent of the spatial coordinates, and indeed we find a time 
derivative of a constraint in one of the linearized field equations (4.23). Also, the differentiated 
constraint (2.17) appears in equation (4.32) in the case that no additional matter other than 
the ponderable aether is present in the universe. 


2.3 Relation to Lorentz-violating dark matter 

Before proceeding with the phenomenology of these theories, it is worth considering what 
physical situations could give rise to a frozen aether theory. In [1], we showed that the frozen 
aether can be thought of either as a fundamental background structure on spacetime, or as 
a type of higher-derivative field theory for four scalar fields, {ip, . In this section, 

we will describe a different situation, where the frozen aether appears as a limiting case of a 
dark matter theory coupled to dynamical Einstein-aether theory. 
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The effects of coupling the aether to dark matter have been explored in [13]. In that 
work, the “pull-back formalism” [18] (see also recent applications in [19, 20]) was employed 
in order to give the dark matter fluid a Lagrangian description. The basic elements of this 
construction are three Euler potentials (f)^, 1 = 1,2, 3, whose gradients 4>l = Vb(/>^ define the 
fluid’s 4-velocity via 


W\' 


r = 


bcda il ij iK 




(t>d ^IJK- 


( 2 . 20 ) 

( 2 . 21 ) 


e 


:0123 _ 


-| In 

Here, |J| = (JVVb) ' , is the anti-symmetric tensor density normalized such that 

= — 1 , and £ijk is a totally antisymmetric set of coefficients with £123 = 1 . 

The action for a perfect fluid respects a symmetry under volume-preserving diffeomor- 

' = 1. The simplest scalar invariant under this 


phisms of the scalars, (j)^ det 




transformation is | J], and the fluid Lagrangian is constructed to be a function of this scalar. 
The form of this function determines the fluid’s equation of state, and it is straightforward 
to see that 


■S'dm = —m 


J d^xv^lJ] 


( 2 . 22 ) 


gives a stress tensor for pressureless dust. 

The lowest order scalar that includes coupling to the aether while maintaining the (j)^- 
diffeo symmetry is 7 = gabU°‘v^. The generalization of S'dm to include Lorentz-violating effects 
simply multiplies the Lagrangian by a function of 7 , 


'S'lvdm = —m 


/ 




(2.23) 


F is further subject to the restriction E(l) = 1 to recover the ordinary dark matter limit 
when the dark matter has no velocity relative to the aether. When 7 — 1 is small, F can be 
organized in an expansion, 

F( 7 ) = l + y( 7 -l) + ... (2.24) 

where Y = F'{1) is a constant. The limit we are interested in is a tight coupling between 
the aether and dark matter, enforced by sending Y —>• 00 . In this limit, and tend to 
align, causing (7 — 1 ) —)• 0 , precluding the need to consider additional terms in the expasion 
(2.24). We will truncate F at order (7 — 1) in what follows. 

Parameterizing the misalignment of by 


ya ^ ^ 


(2.25) 


where u- v± = 0 , we find that since is unit, 7 is related to u_l by 

7 = \/rTKF~l + ^- (2.26) 

Thus for small , 7 will deviate from 1 as the square of the perpendicular velocity |ux|. 

The aether equation of motion receives a correction due to the coupling to the dark 
matter, 


£a - 2\Ua - g-dmYVa = 0 , 


( 2 . 27 ) 






where we have defined /Xdm = m\J\. The u°‘ component determines A, 

A = - ^/idmh^7, (2.28) 

and the z;“ component gives the relation 

= -hdmA^7l^±|. (2.29) 

Since the LHS of this equation should generically be the same order as //dm, this relation 
shows that |u_l| ~ 0(l/y). This also implies that 7 = 1 + 0{1/Y‘^). 

Finally we examine the metric variation: 

T*, + \UaU^ + ^ ((1 - Y)VaVl, + 2yU(„Ub)) (2.30) 

Using (2.28) to eliminate A gives 

+ ^^u'^EcUaUb + ATab, (2.31) 

where AT^b describes the deviation from Einstein-aether theory. The expression for this 
tensor is 


ATaft = ( (7^(1 - ^) + 1^) UaUb 

+ (y (1 - 7 ) + 7 ) 2 u(„u_Lb) + (1 - 5^) v±aV±b) (2.32) 

Now taking the limit Y ^ oo, and using the derived y-dependence of 7 and we hnd 
in the limit the only finite piece is 


AT(jb t (2.33) 

with all other terms falling off as 0(1/y) or faster. We thus recover the frozen aether Einstein 
equation (2.9) with the dark matter energy density /idm playing the role of fi. 

This limit is useful because it allows results about Lorentz-violating dark matter from 
[13] to be applied to the frozen aether theory. In particular, the analytical solutions for cos¬ 
mological perturbations obtained for Lorentz-violating dark matter match on to the solutions 
we consider in this work when the y —)• oo limit is taken. 

3 Small-scale phenomenology 

Here we investigate some astrophysical effects of the frozen aether by considering waves in 
Minkowski space, the Newtonian limit, and spherical symmetry. The most notable effect 
comes from ponderable aethers (those with /i 7 ^ 0) in the Newtonian limit, where /x plays 
the role of a dark matter source. On the other hand, we find that nonzero ^ precludes static 
solutions. The nonexistence of static, spherically symmetric solutions for ponderable aethers 
is a significant challenge, and leads to a puzzle about the hnal state of gravitational collapse. 
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3.1 Linearized waves 

Linearized perturbations around Minkowski space closely resemble the results in Einstein- 
aether theory [21, 22], There are five massless wave solutions, corresponding to spin-2, 
spin-1, and spin-0 excitations. In the background, the aether is aligned with a Minkowski 
time translation killing vector, and /Xa = 0. Hence when perturbations are considered, the 
integration constants must be of the same perturbative order as the metric and aether 
perturbations. 

Using the SVT decomposition of metric and aether perturbations described in Appendix 
B, we find that the scalar variable S obeys the sourced wave equation 

□oS = A7tGn (p + fi + SCaP - 3Ca ^ AtT^^ . (3.1) 

V 2ca + ce J 

Here, Dq = — A is the spin zero wave operator with wave speed 


„2 {^9 + 2Cct)(2 — Ca) 

“ 3ca(l-c<,)(2 + c0)’ 

(3.2) 

S^’^didj is the spatial Laplacian, 


^ _ Go 

9 

(3.3) 


is the renormalized Newton constant for the Einstein-aether Newtonian limit [23], p is the 
matter energy density, and the pressure p and anisotropic stress are defined by decom¬ 
posing the scalar parts of the spatial matter stress tensor by, 

4"^ = P^ij + a) ttl. (3.4) 

Equation (3.1) matches the sourced spin-0 wave equation from reference [22] except for 
the additional term coming from p. The evolution of p given by £uP = 0 becomes at linear 
order in perturbations, 

dtp = 0. (3.5) 

Thus, p \s a. strange type of source that is constant in time. The general solution for these 
wave modes will be the same as in Einstein-aether theory, except the equilibrium value of 
each Fourier mode of E will be shifted by a /i-dependent constant. This zero frequency shift 
would be largely undetectable away from the source. 

The spin-1 analysis produces a similar result. The vector variable V) defined in (B.6a) 


satisfies a wave equation. 


lQ'KG{CaSi + Ki) 

1 l-(l-c^)(l-c^)’ 

(3.6) 

where Di = — A, the spin-1 wave speed is 


„2 1 - (1 - Ca){l - Ccj) 

2Ca(l - c,) 

(3.7) 


Si is the vector part of T^\ and m is the vector part of the source pt. This again agrees 
with reference [22] (see erratum), modulo the extra integration constant k*. As with the 
scalar integration constant, Ki is constant in time, so again only produces a shift in the zero 
point value of V)- 

The spin-2 equations of motion are exactly the same as in Einstein-aether theory, since 
there is no tensor component to pa- 
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3.2 Newtonian limit 


The Newtonian limit is recovered from the linearized equations by requiring slow motion 
of sources so that time derivatives in the linearized equations are suppressed by powers of 
u <C 1. There are two scalar gravitational potentials, given by S and T from (B.9). Equation 
(3.1) then reduces to a Poisson equation 

AS =-47rGAr(/9 +/r), (3.8) 

with Gm defined in (3.3). The 0 component of the constraint (2.11) gives a relation for the 
other gravitational potential, 


A(-2S - CaT) = 87rG(/9 + ^) (3.9) 

which then implies that T satishes the same equation (3.8) as S, so 

S = T. (3.10) 

The Newtonian gauge choice B = U = 0 shows that S and T function as the Newtonian 
potentials (the remaining scalar potential V which functions as a velocity potential for the 
aether is only determined at 0.5 post-Newtonian order). Equation (3.8) demonstrates that // 
will contribute an extra matter component in the Newtonian limit, which could be interpreted 
as the dark matter. In general the distribution of /x need not bear any simple relation to the 
density of ordinary matter p. The value of p, will depend on the details of how gravitational 
collapse occurred from the cosmological evolution. This provides a possible motivation for 
exploring whether p can be an astrophysically viable explanation for dark matter. 

3.3 Static spherical symmetry 

Static solutions for ponderable aethers are somewhat problematic, due a tendency for the 
aether to align with the timelike Killing vector If it does not align, there will generically 
be a flux of p energy density across timelike hypersurfaces generated by flows of This 
will lead to accumulation of p inside of closed spatial regions, which will presumably lead to 
a time-varying source in the Einstein equations. Hence, nonaligned ponderable aethers are 
inconsistent with the existence of a static Killing vector. 

However, when the aether is aligned with p is required to vanish by symmetry. To 
see this, we write and from (2.4) and (2.5) find 

0 = £uPa = N~^£^Pa + Pa£^N~^ - /xValniV. (3.11) 

Since is a Killing vector, £^pa = 0 and £^N~^ = 0. Also, for aligned with the 
acceleration is given by Ua = — ValnN. Since this will not vanish except in the case of 
Minkowski space, (3.11) implies p = 0. 

Imposing the additional restriction of staticity is enough to set the perpendicular source 
density p'^ to zero as well. The alignment of with a static Killing vector implies 6 = 0, 
= 0, and Uab = 0. Also, the relation £uc£ = —u'^a? holds in this case, so examining the 
form of £b given by equation (A.l), we find 

£b oc a?Ub, (3.12) 

so in particular, 

pi = £i = 0. (3.13) 
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This eliminates the possibility of new static, spherically symmetric solutions, outside of those 
that already exist in Einstein-aether theory where /la = 0 [24]. Note that a similar effect 
has been found in projectable Hofava gravity, where accretion of the dark matter integration 
constant prevents static stellar solutions [25]. 

The only possibility left is to drop the assumption that aligns with The general 
argument at the beginning of this section would still apply, meaning that these solutions 
would need fl = 0. However, we no longer can conclude that kinematically vanishes, 
so it may be possible to find fixed threading solutions that generalize the Einstein-aether 
spherically symmetric solutions [24, 26] . This prospect is left to future work. 


4 Cosmological phenomenology 


In this section we look at the cosmological implications of the sources jla- The most interesting 
features come from ponderable aethers, since the background value of will vanish by 
isotropy. The scalar source /i functions as a dark matter component, affecting the background 
FLRW solutions. This is an interesting result, and motivates investigating whether it can 
reproduce other dark matter phenomenology. 

The first issue that must be addressed is determining what fraction of the matter content 
today could be attributed to //. This requires determining initial conditions for //, about 
which, unfortunately, our theory does not have much to say. It seems necessary to simply 
take the initial value of /r as a set of parameters defining the theory. Given a choice of initial 
condition, we can estimate the size of /r today by noting that /r oc o“^. Then the fraction of 
matter comprised of /i today is related to the initial density /Zjnf by 


e 


-3N 



(4.1) 


For an inflationary Hubble rate of 41 ~ 2 x 10^® GeV, a temperature of reheating to be 
Treh ~ 2 X 10^^ GeV, and an initial value for the source to be the same order as the vacuum 
energy density during inflation /iinf ~ pA ~ 2 x 10^*^ (GeV)^, we estimate a value today of 
~ 10“®^. This negligibly small value for effectively reduces the theory to dynamical 
Einstein-aether theory. 

Since the primordial value /i would be washed out by inflation, the only hope for there 
to be observational consequences is for /x to somehow be generated at the end of inflation, 
after which it subsequently satisfies the classical conservation law T„/x = 0. However, this 
situation is fairly implausible as well, since we would need ~ 10“^® to prevent it from 

dominating before matter-radiation equality.^ 

Putting aside the issue of generating the appropriate value of /x primordially, there 
is yet another issue with using /x as a dark matter candidate: it does not lead to structure 
formation. Perturbations to FLRW spacetimes dominated by /x will contain two vector modes 
and one scalar mode, just like the waves in flat space form section 3.1. If /x is the dominant 
form of matter in the universe, it will be solely responsible for the formation of structure 
through gravitational collapse. However, it is already clear that this will be problematic. 
Ordinary cold dark matter, as a pressureless fluid, has growing modes at all scales during 
matter domination. This is not the case for the ponderable aether. Although its stress tensor 

^However, see [6, 27] for proposals for generating fj, at matter-radiation equality in the mimetic theory, and 
[5] for the generation of dark matter in projectable Hofava gravity. 
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behaves like a pressureless fluid in the background spacetime, in general Spu. ~ •So'^/Oee / 0, 
where Sq, dehned in (3.2), is the spin-0 mode speed for Einstein-aether waves. Hence, there 
exists a characteristic Jeans length Aj ~ sq/H (with Hubble rate H), below which the aether 
pressure will prevent gravitational collapse. Since phenomenological constraints coming from 
the absence of gravitational Cerenkov radiation require sq > 1 [28], this means structure 
would not form on any scales smaller than the Hubble length. 

Thus, this model still requires an ordinary dark matter fluid in order for scalar pertur¬ 
bations to grow. The main effect of /r in the presence of ordinary matter is that it changes the 
relationship between % and the background matter energy density p. Thus when is small, 
the cosmological perturbations will be identical to dynamical Einstein-aether theory, which 
have been studied previously [29-32]. Ponderable aethers deviate from Einstein-aether the¬ 
ory only when the universe enters a matter-dominated phase, where may be a significant 
fraction of unity. We will therefore focus on this phase of the universe. The analysis of scalar 
perturbations in this situation reveals growing mode on superhorizon scales, with a growth 
rate depending on the value of This allows us to phenomenologically constrain to be 
less than 0(10“^), precluding its interpretation as a significant source of dark matter. 


4.1 FLRW solutions 


In ordinary Einstein-aether theory, the only effect of the aether on flat ELRW solutions 
is to rescale the gravitational constant appearing in the Eriedmann equations to G^os = 
Go/(l + y) [23]. The ponderable aether gives rise to a new effect; the source p contributes 
additional energy density to the Friedmann equation that behaves like pressureless dust, with 
p oc a~^. 

Working with a flat FLRW metric in conformal time coordinates (B.l), isotropy requires 
that be aligned with the time translation vector dr^, with components 

= a-^6^. (4.2) 

This form is especially convenient for computing Lie derivatives £u and other objects related 
to the aether. For example, the expansion 6 is 


£u\ 9 —4/ —lo 4 I 4o — 

—= a [a d„a + a o„a ) = 3—, 

a 


(4.3) 


with Ti = a'/a, and prime denotes differentiation with respect to p. 

Isotropy imposes that ujab = c^ab = 0, Oa = 0, so that the expansion 6 entirely charac¬ 
terizes the aether. The aether variation (A.l) takes the simple form 


£a = 


(4.4) 


SttGo 3 

the metric variation (A.2) depends only the expansion piece Xab dehned in equation (A.3), 
^ab ~ i6^Go ^^o,b + ^^ab)- The hxed aether Einstein equation (2.9) now reads 

il 

^ - ^>aUb + {2n' + n^)gab] = 87rGo{Tab + pUaUb). (4.5) 

This gives the two Friedmann equations governing the dynamics, 

^- 2^2 ^ ^'kG^cos ^ 

-I- 3p), (4.6b) 
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where p and p are the matter energy density and pressure. 

Hence, the novel feature of this theory compared to Einstein aether theory is the extra 
energy component appearing in the Friedmann equation (4.6a) coming from p. Also the 
requirement that £uP = 0 means the p oc just like pressureless dust. 

In our analysis of perturbations around these backgrounds, we will focus on the matter 
dominated phase where p can represent an appreciable fraction of the energy density of the 
universe. We will use a simple model consisting of the ponderable aether and an additional 
dark matter fluid. Since our focus will be on scalar perturbations, it is useful to assume that 
the dark matter fluid is irrotational. To provide a Lagrangian description of such a fluid, we 
employ the “angel dust” model [33], consisting of a scalar field p and a Lagrange multiplier 
held p with the action 


Sr. = j d^x^g I - l) . (4.7) 

The p equation of motion enforces that Va = Va</? be unit, which also implies that is 
geodesic. The stress tensor is that of pressureless dust, 

T2 = PVaVb- (4.8) 

In homogeneous, isotropic cosmology, and since p for the dark matter vanishes, 

the solution to the Friedmann equations (4.6) is a matter-dominated universe with a (x rf, 
% = 2 / 77 . The dark matter helds satisfy 

p' = a, (4.9a) 

with p. Note that during matter domination, remains constant. 

4.2 Scalar perturbations 

Now we turn to deriving equations for linearized perturbations about this solution. The 
metric and aether perturbations are parameterized using the SVT decomposition described 
in Appendix B. The equations of motion can be derived by expanding the actions (2.1) and 
(4.7) about the background solution to quadratic order in the perturbative variables. Varying 
the resulting quadratic action yields the linearized equations of motion. 

For the ponderable aether, we must also ensure that the linearized aether equation of 
motion is not imposed by this procedure. These equations result from varying the aether 
perturbations T, U and [/*. To avoid enforcing these equations of motion, we can gauge-hx 
these quantities at the level of the action, so that their variations are not included. Choosing 
T = U = W = 0 simply selects the gauge where = a~^6Q. This also allows us to solve the 
unit constraint for the goo metric perturbation, and gives A = 0. When this is substituted 
back into the action, the A term ( 2 . 2 ) drops out. 

These conditions in fact do not fully hx the gauge; there are still coordinate transforma¬ 
tions that leave this form of u" invariant. These correspond to shifts of the time coordinate 
along each of the aether threads, 1 1 —>• t + f{x^), and thread-preserving reparameterizations of 
the spatial coordinates, x* x^{x^). Rather than trying to fully eliminate this residual gauge 
ambiguity, we can instead focus on gauge invariant combinations of the remaining nonzero 
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variables. These combinations are given in equations (B.5), and in our choice of gauge take 
the form 


Y. = -C-'HB, (4.10a) 

+ (4.10b) 

V = E', (4.10c) 

and = denotes equality only in this gauge. 

The matter field ip and p also receive perturbations, 

p = p + 6p (4.11a) 

(f = ip + 6(p, (4.11b) 


where the bar denotes the background value of the fields, given by (4.9). The useful gauge- 
invariant quantities associated with these are 


6ip 


ip> 


M = n[B + u + 

Z = E + M = -C + n 




= n{B + 

Sip 


Sip 

ip' 


ip' ’ 


(4.12a) 

(4.12b) 

(4.12c) 


The variable Z has the geometric interpretation of the curvature perturbation on surfaces of 
constant ip. When we consider only scalar perturbations, the aether is orthogonal to surfaces 
of constant B + U, so that S gives the curvature perturbation on these surfaces. The variable 
M is then related to the difference in conformal time between these surfaces [31]. The final 
variable is simply the matter density contrast in the flat-slicing gauge, where (7 = 0. 

It is convenient to work with the spatial Fourier transform of the perturbative variables, 
defined by 


X{-K,r]) = 


d^k 

(27r)3/2' 




XUv)- 


(4.13) 


From now on, we will work with the Fourier-transformed variables X]^ and drop their explicit 
dependence on k. 

The aether action (2.1) is now expanded to quadratic order. The details of this cal¬ 
culation for the aether quantities are given in Appendix C. Using the background equation 
2'H' -l- B? = 0 and integrating by parts, the action takes the form 


oae _ 

IQttGq 


where 


kl 


dr] 


(2 + ce)^,2 _ 2(1— c4.{4E'f - Cak‘^a-^{aBy^ - 2k^{Ty - 2ST) 


12 


h = 6C- 2k^E. 


(4.14) 

(4.15) 
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The matter action (4.7) at quadratic order is 


cm _ 

*( 2 ) - 


IdvrGo 


— m) 


^Px.J M 


^ V - k 
L P 




(V 


r\2 


{if' 


n ''\2 


+ k^ 


6f^ 

w? 


(4.16) 


with 

m =-y+ (l + y) . (4.17) 

Here we have used equation (4.9) to express p and (p' in terms of a and %. 

The total quadratic scalar action S^ 2 ) + 'S'( 2 ) is thus a function of five independent 
variables {h, E, B,6p,6if), since S, T, and M are simply functions of these. We expect there 
to be two dynamical modes for this system corresponding to the aether scalar mode and 
the angel dust mode. The equation from varying p takes the form of a constraint, and from 
the discussion in section 2.2 there are two other constraints with undetermined integration 
constants, given by equation (2.18) and the scalar part of equation (2.19). These will reduce 
the original five variables down to the two dynamical ones. 

The 6p variation gives the first constraint,^ 

Sip' = 0. (4.18) 


This helps simplify the equations; for example, the hSip' term in (4.16) no longer contributes 
to the h equation. Varying now with respect to Sp yields another constraint, which in terms 
of gauge invariants reads 

“ f) “ = "■ 

Since does not enter any of the other equations, this equation simply determines it once 
the solutions for M and V are known. 

The remaining parts of the action now only depend on three variables h, E and B. It 
is more convenient to change variables from h and E to 

C = \eE + \h, (4.20) 

6 6 

K = \eE - \h. (4.21) 

6 6 

After doing this, the X-dependent terms of the action take the form (from now on we ignore 
the overall l/(167rGo) prefactor, which is common to all terms) 

= j dp [{ce + 2c^)K'^ - 2(4 + ce - 2 c„)C'K') , (4.22) 

Here we see that K is cyclic, so its equation of motion will be a total time derivative. This 
is in fact the constraint expected from equation (2.19), which after integrating contains an 
undetermined constant ps corresponding to a source density, 

{ce + 2c^)K' - (4 + ce - 2c,)C' = a-^.- (4.23) 

^The gauge-invariant version of this equation is ^- A = 0, but T = 0 for our gauge choice. 
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This can be cast into a gauge-invariant form, using K' = ^k'^V — C', and C = —H — |'HC, 
which gives 

— {cg + 2ccr)k‘^V + 2(2 -|- cg)^ = (a — 3(2 + cg)T-LC^ . (4.24) 

O 

The gauge-invariance of the right hand side follows from viewing as arising from the 
scalar part of the spatial component of the constraint violation, = di^g- Under a gauge 
transformation generated by this changes by where the background 

value of is determined via the background Einstein equation to be 3(2-|-C0)a^'H^. Gauge 
invariance then follows by noting that C transforms as 6C = Equation (4.24) determines 
V once a solution Z is known, since E = — ^ ^ U determined by Z up to a term 

proportional to 6ip, which by (4.18) is just an integration constant. 

Equation (4.23) can be used to eliminate K' from the action (4.14). When this is done, 
the fig piece decouples, and up to total derivatives the action for the two dynamical modes, 
including the piece from (4.16), can be written in terms of gauge invariants. 


Qdyn _ 
*^(2) “ 




-|- a^k 


2l2 


2^2 -P 4—Z' -h 3mM 
T~L 



where 


6(2 + cg)(l-0 

Cg 2Cfj 


(4.25) 


(4.26) 


Erom this action we can derive the two dynamical equations for Z and M. Using the back¬ 
ground equation a = rf, they take the form 


pZ + —Z + yM -2Z + 5M = 0, 

y 


2^m + ^m--z + 

2 y y 


3{ca - 2m) 


M = 0, 


(4.27a) 

(4.27b) 


where y = krj and dots represent differentiation with respect to y. 

These equations describe how the two scalar modes evolve, and the dependence on 
appears through the parameter m in (4.27b). An important set of solutions to these equations 
comes from the long wavelength limit y ^ 1, or equivalently fcphys = k/a ^ H = 'H/a, where 
equation (4.27a) reduces to 

/3(^Z + Jz^=0, (4.28) 

and equation (4.27b) remains the same. These equations now admit four independent power 
law solutions. The first two solve (4.28) with 


Zocy''±, g± = 0,-3. 


(4.29) 


For the = 0 solution, M is simply given by M = 0. For q- = —3, M obeys the same 
power law as Z, with an amplitude determined by equation (4.27b) to be 


M_ _ 1 
Z_ m 


(4.30) 
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Then using equations (4.23) and (4.19), we find that the matter density contrast 6^ evolves 
as 6^ (X so that the growing mode satisfies oc a. This growth rate is the same as 

cold dark matter during matter-domination. Thus, we see that the matter decouples from 
the scalar aether perturbations on superhorizon scales. Once the mode enters the horizon, it 
will generically mix with the aether perturbations and lead to decaying oscillatory solutions. 

The other two solutions set Z = Z = 0, which is preserved by (4.28). Then (4.27b) has 
power law solutions M oc with 


P± 


-l±Wl + 


24 0^(2 -h Cg) -Cg-Ca 


25 


(4.31a) 


5 

2 


-1± Jl + 


48Go / I 1 - 


25ca \Gn 


Gr 


(4.31b) 


Results that reduce to this exponent when = 0 have previously been derived in Einstein- 
aether theory [31], generalized Einstein-aether theory [34], and Hofava gravity [35]. Eor 
nonzero values of 0^, these exponents can also be derived from the analysis of Lorentz- 
violating dark matter [13], by taking the parameter Y —)• oo, as described in section 2.3.^ 
Another interesting set of solutions occurs when the universe is dominated by the pon¬ 
derable aether, i.e. 1, since in this case the order of the equations reduces to a single 

dynamical scalar. To see this, we add (4.27a) to times (4.27b) and multiply by y, resulting 
in 


A. 

dy 


^^{y^Z) - y^Z + ^^{y^M) + y^M 


dy 


4 dy 


3(m — l)yM = 0, 


(4.32) 


which can be immediately integrated when m = 1. This gives a constraint equation that 
can be written entirely in terms of T, and T, and the integration constant resulting from this 
procedure is simply the perturbation to the source y. 

Solving the constraint for T, and inserting the result into (4.27a) results in a single 
evolution equation written entirely in therms of S, 




\. 2u 1 


4 + 2 1 

yS -b 

[ y^ + v] 



2 2 20 , 8(3/3+ 1^) 

Ca Ca(?/2 + u) 


s = 


-2 4 

/3 y"^ + u 


c, 


(4.33) 


4/3 

where v = — and Sq is the spin-0 aether wave speed from (3.2). 

Ca 

This equation does not have an exact analytic solution; however, we can describe its 
behavior in limiting regimes. At short wavelengths, s^y ^ 1, the equation takes the form of a 
Bessel equation with a time-dependent, decaying source. Solutions thus behave like damped 


^ Note that the exponents in (4.31) diverge in the limit that Ca —>■ 0, and, as explained in section 4.3, 
small values of Ca are phenomenologically preferred. On the other hand, if we set Ca = 0 in the action, the 
M field equation (4.27b) becomes a constraint relating M to Z, eliminating one of the scalar modes. Hence, 
we see the Ca —>■ 0 limit of (4.31) does not match the Ca = 0 solutions. This is related to the strong-coupling 
issues that arise in Hofava gravity and Einstein-aether theory in the limit that the parameters Ci become 
small [36-38]. In particular, the effective field theory becomes strongly coupled at energy scales of the order 
A = yTTAfp. For Hofava gravity, the higher spatial derivative terms can provide a UV completion of the IR 
theory defined by the action (2.1) [37]. It is unknown whether a UV completion exists for Einstein-aether 
theory. Nevertheless, the strong coupling scale A can remain close to the Planck scale even for relatively small 
values of Ca. 
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oscillators with a shifted equilibrium value. In the opposite regime squ <C 1, we obtain power 
law solutions describing growth on superhorizon scales. The exponents are given by equation 
(4.31) after setting 17^ = 1. 

4.3 Constraints on 

The superhorizon growth exponent from (4.31) allows us to constrain the size of 17^. 
The quickest constraint comes from requiring that the potential M remain linear throughout 
the epoch of matter domination. This requirement is justified since M contributes to the 
gravitational potential; it is related to the usual Bardeen potentials <h and T [39, 40] by 


^ = z -M + nv, 


4* = a ^ 



(4.34a) 

(4.34b) 


Since these gravitational potentials are linear except on scales as small as stars or black holes, 
we require that M should be linear as well. Although this leads to a rather modest constraint 
on p_|_, it translates to a much more stringent constraint on 17^, allowing us to conclude that 
/i cannot comprise a sizable fraction of the dark matter. 

When the post-Newtonian constraints are satisfied by the parameter choice (4.37), the 
analysis of primordial perturbations for Einstein-aether theory [31] reveals that both super¬ 
horizon modes found in section 4.2 should have comparable amplitudes of the order of 10“^ 
at the time of radiation-matter equality. The superhorizon growing modes would have an 
amplitude today of 

Mo = Meq(ao/aeq)'’+/'. (4.35) 

With a redshift at equality of Zeq ~ 3500 and initial amplitude Meq ~ 10“^, Mq will remain 
less than unity if 

P+ < 3. (4.36) 

To translate this to a bound on 17^, note that solar system tests of post-Newtonian 
parameters require Gn and Geos to be very nearly equal [41, 42]. By examining the form of 
in (4.31b), we see the Geos and Gw terms cancel when the post-Newtonian constraints 
are satisfied, and the remaining term is proportional to Q^/ca- Since Ca is phenomenologi¬ 
cally constrained to be small, 17^ must be similarly small to produce an order unity growth 
exponent. 

More quantitatively, the preferred frame post-Newtonian parameters ai and ol 2 can be 
set to zero in Einstein-aether theory by choosing [41] 


-=- = U-+-) 

Ca ^ \Cf7 C(jj J 


(4.37) 


The Cerenkov constraint imposes that both and are positive, and observations of binary 
pulsar decay rates further imposes that [43] 

< 3 X 10-2, 

Co; < 3 X 10“^, 

which, combined with (4.37) gives 

Ca < 5 X 10-3. (4.39) 


(4.38a) 

(4.38b) 
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Now, given a bound on the exponent < b, and using the above restrictions on the 
aether coupling parameters, this translates to a bound on 




< 

rs_/ 


6(6 + 5) 
12 


(4.40) 


Taking 6 ~ 3, this leads to 

% < 10-^ (4.41) 

Such a small value of 11^ precludes it making any appreciable contribution to the dark matter 
in the universe. 

The bounds on and could be improved by considering the contribution of the 
growing modes to the integrated Sachs-Wolfe effect in the CMB. The time varying gravi¬ 
tational potential 4> — T during matter domination coming from the growing modes causes 
changes in photon wavelengths as they free stream from the surface of last scattering. If the 
growth rate is too large, this will lead to an excess of power in the low I multipoles of the 
CMB power spectrum. 

Finally, we note that this constraint only applies for Ca finite and nonzero. In particular, 
it does not apply to projectable Hofava gravity or the mimetic dark matter theory. The 
growth exponent for those theories is obtained from the ponderable aether result by taking 
the limit Cq —)• oo in (4.31). This results in p_|_ = —2 and = —3, so that the isocurvature 
modes are decaying on superhorizon scales. These decaying modes would have little impact 
on CMB and large scale structure observations. 


5 Discussion 

This work has explored several consequences of treating the aether vector u°‘ in a Lorentz- 
violating theory as nondynamical. The resulting ponderable aether theory differs from dy¬ 
namical Einstein-aether theory by the presence of source densities in the aether equation of 
motion. The source densities take the form of a covector density that is conserved along 
flows of They also have an interpretation in terms of the violated initial value constraints 
of Einstein-aether theory (section 2.2). 

The most striking phenomenological feature of the source densities for ponderable 
aethers is that the scalar density /2 can mimic the effects of cold dark matter in FLRW 
cosmologies (section 4.1). However, beyond homogeneous, isotropic cosmolgy, \i differs from 
cold dark matter on several fronts. First, as argued in section 4, the significant pressure 
in the aether fluid prevents gravitational collapse on scales shorter than the aether sound 
horizon, which is larger than the Hubble radius. Thus it does not lead to the structure 
formation observed in our universe. Furthermore, perturbations to the FLRW backgrounds, 
considered in section 4.2, show that the ponderable aether produces a problematic growing 
mode on super-Hubble scales unless is relatively small. This precludes its interpretation 
as a significant component of the dark matter. On smaller scales, /r does not allow for static, 
spherically symmetric solutions, so that it is unclear what the hnal state of gravitational 
collapse would be. 

On the other hand, considerations of cosmological inflation provide a different perspec¬ 
tive on this theory. Under reasonable assumptions for the value of /i at the beginning of 
inflation, the requisite 60 e-foldings is more than enough to dilute it to a negligible value 
by the present era. This is similar to the inflationary solution to the monopole problem 
[44, 45], where in our case is playing the role of an exactly stable massive particle. With 
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H diluted to nearly zero, the theory essentially reduces to dynamical Einstien-aether theory. 
Thus, a universe described by Einstein-aether theory could equally well be interpreted as a 
universe with a nondynamical, Lorentz-violating background structure, whose effects have 
been washed away by inflation. In this case, the line between the diffeomorphism-invariant 
and noninvariant theories is not sharp. 

Ponderable aethers are a vector field generalization of the mimetic dark matter con¬ 
struction [6], related to the model proposed in Ref. [12]. By taking the aether field to be 
the gradient of a scalar (as was considered in [46]), we obtain the IR limit of projectable 
Hofava gravity [8, 9]. The aether kinetic terms then become the additional gradient terms 
that appear in extensions of the original mimetic dark matter theory [27, 47]. Note that 
these kinetic terms would be expected from an effective field theory point of view, and are 
also necessary when the aether is nondynamical to avoid overly constrained gravitational 
dynamics [1]. We note that the constraint derived for ponderable aethers in section 4.3 does 
not apply for the mimetic and projectable Hofava theories, since these theories come from 
a Ca —)• oo limit of the ponderable aether theory, where there are no growing superhorizon 
modes. 

Although the ponderable aether is highly constrained phenomenologically, it never¬ 
theless presents some directions of theoretical interest for future work. The cosmological 
constraints derived in section 4.3 were rough estimates, intended to only bound enough 
to exclude its interpretation as the dark matter. More detailed analysis of the cosmological 
implications could be performed, and likely could improve these bounds on significantly. 

The non-existence of static, spherically symmetric stars in this theory is a distinctive 
feature, and presents a puzzle regarding the final state of gravitational collapse. This could be 
resolved by examining numerical simulations of gravitational collapse, along the lines of [48] . 
Another possibility would be to look instead for stationary, axisymmetric solutions, which 
allow [I to orbit a central star or black hole without accreting. Another line of investigation 
could focus on fixed threading theories, where /r = 0, but is nonzero.^ In that case, by 
choosing the aether vector to not align with the Killing vector, there may be solutions that 
generalize Einstein-aether black holes [26] with sourced spatial aether equations. 
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'‘Note that this is just a special case of the fixed aether theory where the initial value of ^ is set to zero 
by hand. One can justify such a choice if the background structure is a fixed threading, as opposed to a fixed 
aether [1]. 
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A Formulae for equations of motion 


1 SS 

The variational derivative of (2.1a) with respect to the aether, Sf, = —p=-j—is given by 




£h = 


1 


ce 


+ Cai-o^Ub + 2a^u}cb + 2V{uicab])) 


SttGo L 3 

+ Caji'V‘'u}cb - a'^UJcb) + Cct(V^CTc6 + O^cb) 


(A.l) 


The metric variation of (2.1a) with held hxed, = 


1 6S, 


CIS _ 

^ab — 


-1 


(Gab + + CaAab + C (jSab T Ca]'Vab 


takes the form 


IGvrGo 

where Gab is the Einstein tensor, and the other tensors appearing are 
A’afe = gab (^2^^ 

Aab —Qj gab T T 2Vc(m (^(b CL U(^b')'^a) 

Sab — 2 *^ 9ab ^bc T ~^(^CTab T £uCCab 

Fab — 2^“^^ ^bc T 2(V Cl (^c{a)'^b) 


(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 


B Perturbative variables 

In both the analysis of linearized perturbations to Minkowski space considered in section 
3.1 and the analysis of cosmological perturbations in section 4.2, it is convenient to use a 
scalar-vector-tensor (SVT) decomposition of the perturbative variables. The reason is that 
the modes of different spins decouple in the linearized theory, and the SVT decomposition 
makes this explicit and also allows for the construction of gauge invariant variables. Here, 
we describe the decomposition used throughout this paper, following closely the conventions 
of [40]. 

Starting with cosmological perturbations, the flat background metric with conformal 


time is 

ds'^ = a^{g){dg^ — Sijdx^dx^), 

(B.l) 

and isotropy requires 

a background value for the aether to be 



= (o"\0,0,0). 

(B.2) 

The perturbations to 

the metric are parameterized by 



500 = a^(l + 2A) 

(B.3a) 


gQi = -a^idiB + Bi) 

(B.3b) 


gij — CL {dij -|- hij') 

(B.3c) 


hij = 2G5ij -\- 2didjE 2d(^{Ej'^ 2Eij, 

(B.3d) 
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and the aether perturbations by 


(B.4a) 

(B.4b) 


= a-^{5^^djU + U^). 

As is standard in the usual SVT decomposition, perturbed quantities carrying an index are 
divergenceless, i.e. diB^ = diE^^ = 0, and the tensor perturbation is also traceless, 6^^Eij = 0. 
Note also that we use 6ij to raise and lower indices of perturbed quantities. 

From these we can construct gauge-invariant variables that appear in the linearized 
equations of motion. The following gauge-invariant scalars are useful, 


E = -C -n{B + U) (B.5a) 

T^A + iB + uy + niB + U) (B.5b) 

V = E' + U (B.5c) 

= -^((|) +5] + t), (B.5e) 


where ' denotes differentiation with respect to rj, and T-L = a'/a. A set of gauge-invariant 
vector variables are 


Vi — Bi -\- Ui (B.6a) 

Wi = E' + Ui. (B.6b) 

Also note that the tensor Eij is gauge-invariant. 

The SVT decomposition of perturbations to Minkowski space can be obtained from 
these formulae by setting a = 1 and = 0. In particular, the metric (B.l) becomes 

ds^ = dt^ — 5ijdx^dx\ (B.7) 

where we have renamed ?? —)• t, and the background aether is just 

u" = (1,0,0,0). (B.8) 

For reference, we also list the relevant gauge invariant scalars: 

E = -C (B.9a) 

T = A + B' + U' (B.9b) 

V = E' + U. (B.9c) 

The expressions for the vector gauge-invariants (B.6) and the tensor gauge-invariant are 
unchanged. 

C Perturbations in the aether gauge 

The aether gauge choice greatly simplifies the computation of the quantities appearing in the 
aether action (2.1) in perturbation theory. In this section, we describe the calculations that 
lead to the quadratic action for cosmological perturbation given in (4.14). 
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The scalar part of the perturbed metric determinant is 

V^{s) = (^ + -jh + , (C.l) 

where h = 5^^hij = 6(7 + 2AE, h 2 = and A = 5^^didj is the spatial 

Laplacian. 

This leads to the expansion 

(C.2) 

which gives the quadratic order cg piece of (2.1), 

iv^o^) (2) = + 9n^h2 + mh’^ , (C.3) 

The last term can be integrated by parts in the action to give 

+ 2n')h^ , (C.4) 

and so in the special case of a matter-dominated background that we are considering, the 
second term drops out. 

Since the other parts of the Lagrangian (2.1) vanish in the background, the acceleration 
and shear only need to be computed to linear order, and the twist contains no scalar piece. 
Also these quantities are orthogonal to which in the aether gauge means they have no 
0-component. The acceleration is 

Oq, = £uUa = —5\a~^{adiB)'. (C.5) 

The shear takes a simple form by first defining the spatial metric 


and then using the equation 

Contracting these quantities using the background metric yields the and Ca pieces of 
(4.14). 
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